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Abstract 



| We refute past suggestions that the data on soft diffraction dissociation demand modifications to 
■ conventional Regge theory. 
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Hard diffractive processes have been an active area of study at HERA and the Tevatron, and will 
CN \ continue to be so at the CERN LHC, but one must be uneasy about all theoretical analyses if we do not 
even understand soft diffraction. The dataW PI from the CERN ISR for soft diffraction dissociation, 

^> : pp^px (i) 

CS ' with the proton emerging with almost no energy loss, have been shown PIM to be in good agree- 
ment with the triple- Regge description^]^]. However, the higher-energy data from the CERN SPS 



collider [^J [8] do not show the increase of the cross section with energy expected from simple fits^l. 

This has led certain authors t 5 l t 9 l to introduce unconventional features into the data analysis, which 

yield energy dependence in conflict with the standard notions of Regge theory^] [6] _ T n this paper, 

we show that this is unnecessary and that conventional Regge theory gives an adequate description of 

Q" 1 - the data. 
<D . 

Let £ = 1 — xp be the fractional energy loss of the emerging fast proton (or antiproton) and t the 
squared momentum transfer to it. According to standard Regge theory, when £ is small enough and 
the energy large enough the reaction (1) is generated through the very fast proton "radiating" a 
5— i ' pomeron, which hits the other proton and breaks it up. This gives 

a | 

^ = D («)^(M=, t ) (2) 

where M is the invariant mass of the system X in (1), 

M 2 ~ £s (3) 

This factorisation of the differential cross section into a product of a "pomeron flux" factor D(t, £) 

and a pomeron-proton "cross section" a Fv is ambiguous; different definitions in the literature W 
differ by a factor of ^7r in the definition of D(t,£), with a corresponding difference in the definition of 

a IPp . Our own definition isW 

D(t,0 = ^(F 1 (t))*e- 2aF(t) (4) 
* email addresses: ad@a35.ph.man.ac.uk, pvl@damtp.cam.ac.uk 
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Figure 1: triple- reggcon diagram 

where is the proton's Dirac form factor, 

aip(t) = 1 + ep + a'jpt 

ejp = 0.0808 a'jp = 0.25 GeV -2 (5) 

is the pomeron trajectory and (3p « 3.5 GeV -2 is defined such that the pomeron-exchange contribution 
to the total cross section is^ 

a pp (s) ~ 18 /^(c^s)^ cos(±7r eiP ) = 21.7s - 0808 (6) 

in mb-GeV units. 

Experimentalists commonly integrate d 2 a/dtd£, = D(t,£) a 1Pp (M 2 ,t) over all t and down to some 
value M MIN of M, and refer to a total diffractive cross section 

/*£=£max 

a DIPF ( S ;M MIN ,W) = / dtdid 2 a/dtdi (7) 

J M=M M in 

The contribution (2) makes <r DIFF (s; M MIN ) rise with s almost as fast as s 2ep . As this is faster than 
the total cross section, it has been claimed^] that shadowing effects are important. However, at 
present energies <r DIPF (s; M MIN ) is much less than the total cross section, and so it seems to us likely 
that shadowing is not yet an issue. This is controversial: for a number of reasons we believe® that 
shadowing is still unimportant also in forward elastic scattering and in the total cross section, but this 
view is not universal. 

Being a hadronic cross section, a ]Pp (M 2 ,t) is expected to be dominated by pomeron exchange when 
the corresponding centre-of-mass energy M is sufficiently large, and then it becomes proportional to 
(M 2 ) ejp , in analogy with the behavior s eF of a pp (s) given in (6). When the centre-of-mass energy 
t/s for a pp (s) or M for a Fp (M 2 ,t) is not large enough, pomeron exchange alone is not sufficient, 
and nonleading Regge exchanges become important. The dominant such exchange in a pp {s) is ji 
exchange®, though p, u and a2 exchanges are present too. Each of these contributes to a pp (s) a term 
behaving approximately as 1/y/s. When this is added to the slowly-rising pomeron-exchange term 
(6), the result is that, as s increases, a pp {s) initially falls with increasing s, and then begins to rise. 
We expect that /2 exchange will similarly contribute to a ]Pp (M 2 ,t), so that it too falls initially with 
increasing M, and then begins to rise. Note also the very important point that, for s less than some 
value sq, these exchanges are not sufficient. For most hadron-hadron total cross sections the value of 
i/s ^ s found to be about 4 to 6 GeV. Similarly, we should expect the simple exchange theory to 
break down for a 1Pp (M 2 ,t) when M is less than some value Mq. One cannot know what the value 
of Mq is; we will be fairly bold and assume a value 3 GeV. This limitation is important when one 
is interested in the total diffractive cross section cr DIFF (s; M Mm ). Usually the value of M Mm that is 
chosen is less than the likely value of Mq, in which case the integration over M extends into a region 
where there is poor theoretical understanding. 
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According to what we have said, when M is large enough for a Fp (M 2 ,t) to be approximated by 
pomeron and other reggeon exchanges, d 2 a/dtd^ in (2) is the sum of triple-reggeon contributions. 
See figure 1. In this figure, the two upper zigzag lines are the pomeron "radiated" by the very-fast 
proton. The lower one is the pomeron or reggeon exchanged in a 1Pp (M 2 ,t). When £ is not very small, 
there are other contributions from a reggeon other than the pomeron being radiated by the very-fast 
proton. Also, these additional contributions can interfere with the pomeron-radiation term. So we 
obtain a sum of terms of the type shown in figure 1: 

pp pp hh hh hP JPh hP "P 

P h P h P P h w U 

The last four terms are examples of interference terms. All the terms but the last contribute equally 
to pp — > pX and pp — > pX or pp — > pX. However the present data turn out not to be sensitive 
to these terms so we treat pp — > pX and pp — > equally There are also contributions from 
and p exchange though, as we have said, these couple more weakly to the proton or antiproton. For 
simplicity we assume that the reggeon exchanges are degenerate. At small t, it is important also to 
include 

7T7T 7T7T . . 

JP A (9) 

Altogether, there are a very large number of terms. For each there is a triple-reggeon vertex, the 
blob in the centre of figure 1, whose dependence on t is completely unknown. With so much freedom, 
it is surely possible to obtain a good description of the data without having to introduce new and 
unconventional devices into the fit. 

There is a huge amount of data on soft diffraction dissociation. If one plots them one finds that a lot 
of them are less than good. ISR and SPS collider data have resolution problems for small £ and so we 
have restricted our use of these data to £ > 0.02. It is a great pity that some of the SPS collider data 
for d 2 a/dtd( t survives only in integrated form. It is unfortunate that none has been published from the 
Tevatron, though Goulianos and Montanhat 1 ^ have reconstructed CDF data at —t = 0.05 GeV 2 using 
a particular model. There are some apparently very accurate fixed-target datat 16 ! though this data 
set has certain internal peculiarities and at some values of t its normalisation disagrees with that of 
references [1] and [2], which agree with each other. Figure 2 shows the fixed-target data at two values 
of s and illustrates that any successful fit must have quite a complicated structure. The differential 
cross section at fixed £ decreases with increasing s in this energy range, which is the opposite behaviour 
from that of the term in (8). At very small t the differential cross section is very steep and the 
data cannot be parametrised with a single simple exponential e , although this fits satisfactorily at 
larger t. Simple exponential fits are often used by experimentalists to extrapolate their data beyond 
the region where they are measured, but this is known to be risky. For example, it is knowiJ^ that 
the slope b decreases rapidly as £ increases. 

We have chosen to use collider data from reference [2] at s = 2880 GeV 2 . They are relatively smooth, 
have quite small error bars, and extend over a wide range of t, from — t = 0.11 to 0.84 GeV 2 . Note, 
though, that it may be unsafe to extend the simple Regge-exchange picture to such a large value of t. 
As we said, to avoid resolution problems we restrict ourselves to £ > 0.02 with these data, and at this 
high energy this corresponds to M > 7.6 GeV, comfortably satisfying our condition M > Mq when 
M = 3 GeV. To give leverage in s, we use SPS collider data from reference [7] at —t = 0.55 and 0.75 
GeV 2 . We use also all the fixed-target data from reference [16]; they are at s = 262, 309, 366, 565, 741 
GeV 2 and small t, -t = 0.036,0.075,0.131 and 0.197 GeV 2 . For these low-energy data the constraint 
M > Mq is important. We use also the reconstructed CDF datat 15 ! at -t = 0.05 GeV 2 . 

A simple least-x 2 fit to these data is not adequate because the error bars on the ISR data are so much 
larger than on the fixed-target data. So we artificially reduce the error bars on the smallest-t ISR data 
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when we make the fit. Also, if we use only the reconstructed CDF data at the one value of t for which 
it is available, when we integrate the fit over £ its t dependence is in conflict with what we believe 
to have been measured by CDF. So we have further reconstructed some CDF data for ourselves, at 
-t = 0.01 and 0.09 GeV 2 with ^ = 1800 GeV. 

Our purpose is to show a reasonable fit to the data within conventional triple- Regge theory, to demon- 
strate that it is not necessary to introduce any new and unconventional ideas. In order to keep things 
as simple as possible, we have just used the first seven terms in (8), together with the pion terms (9). 
The latter contribute only at very small t and are known, essentially without any free parameters: 

^ = D*(t,0<r* p (M*,t) (10) 

with 

n^a ft = 

16vr 2 (t - ml) 2 



^-13.3 a n (t) = a' 7T (t-ml) (11) 



An 

For the trajectory slope we use a' n = 0.93 GeV -2 . For the form factor G(t), which takes account 
of the proton wave function, we use the Dirac form factor; this is surely not correct, but it will not 
matter because d 2 a v /dtd£ is significant only at very small t. For the same reason, we ignore the t 
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Figure for (\jii)d 2 o j dt d£ in mb GeV 2 plotted against £ at various energies and values of t, 

compared with the fit. In each case, only data points to the right of the vertical line are used in the fit. 
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Figure 4 continued: open points are from reference [1], black points from UA4[ 7 L In each case the upper curve 
is at the higher energy. 
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Figure 5: (a) reconstructed CDF data^J for (l/ir)d 2 a/dtd^ in mb GeV -2 plotted against £ with our fit; (b) 
the rise of (l/ir)d 2 a/dtd^ with y/s at —t = 0.05 GeV 2 (upper pair of curves) and 0.5 GeV 2 . In each case the 
upper curve is for the higher energy. 



dependence of a np (M 2 ,t) and use its on-shell formal 

a^{M 2 ,t) = 13.63 (M 2 ) - 08 + 31.79 (M 2 ) -0 - 45 

We have used a linear f% trajectory^] 

e R = -0.45 a' R = 0.93 GeV" 2 

and an f 2 flux factor 



(12) 



n2 el-2a R (t) 



D H (t^) = ^(F 1 (t)) 2 i 

with F\(t) again the proton's Dirac form factor. We also use 

a Fp (M 2 ,t) = VjpjpjpH) (M 2 ) om + Vjpjp R (t) (M 2 )" - 45 

a R v(M 2 ,t) = V RR]P (t) (M 2 ) - 08 + V RRR (t) (M 2 )- - 45 
where the Vi(t) are the triple-reggeon vertices. 



(13) 
(14) 



(15) 
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Figure 6: (a) da/dt in mb GeV -2 plotted against \t\ for various values of y/s; (b) <r DIFF plotted against y/s. 
Both are for £ integrated up to 0.05 and down to M 2 = 1.4 GeV 2 . In (b) the middle curve is for —t integrated 
down to 0.036 GeV 2 , the lower curve is the contribution to it from F p and the upper curve is an estimate for 
—t integrated down to 0. The data are from the table in reference [15]. 

It is well testedt 18 ! in pp elastic scattering that it is the Dirac form factor that takes account of 

the proton wave function in the coupling of the pomeron to the proton. This is why we have included 
Fi(t) in the pomeron flux factor D(t,£) in (4). The coupling of the reggeon to the proton is less 
well determined. Nevertheless, we are free to include the same form factor in the reggeon flux factor 
Dn(t,£) in (14), because in the differential cross section it is multiplied by triple-reggeon vertices. 
There is no knowledge of how these four triple-reggeon vertices depend on t, and the only way we can 
determine them is to fit their products with the flux factors D(t,£) or D R (t,^) to experiment. To 
this extent, what we choose for the t dependence of the flux factors is just a matter of convention; it 
is only when one introduces some assumption about the t dependence of the triple-reggeon vertices, 
or of a 1Pp (M 2 ,t) and a Rp (M 2 ,t), that the definitions of the flux factors take on some independent 
significance. 

Previous fits™" [19] nave included large interference terms, and so we have included such terms 
maximally. By this we mean 




The cosines arise from the Regge signature factors®. 

We have experimented with various simple forms for the triple-reggeon vertices and have arrived at 
VAt) = Ai (— *—) Xr i = P1P1P,1P1PR,RR1P,RRR (17) 

These are multiplied by the flux factors, which include the couplings 0jp or (3 2 R . We have varied the 
12 parameters f3 2 Ai, /if and Ai to make our fit. It has Xrrr quite large and negative, and Xppr also 
negative. One should not conclude from this that the data require that the corresponding vertices 
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diverge at t = 0, because the smallest value of t for which there is data is —t = 0.036 GeV 2 and it is 
only for \t\ larger than this that our fit has any meaning. As may be seen in figure 3, our vertices all 
behave sensibly. The reason that Vrrr needs to vary fairly rapidly at small t may be seen in figure 2: 
at least at low energy, the data fall steeply with \t\ at very small t and then flatten off. It has been 
known for a long timet 2 ^] that Vpppft) is small at small t. 

In figures 4 and 5a we compare our fit with the data we have used, except that we have not shown 
the s = 741 GeV 2 data from reference [16] because they are rather ragged. Note that our fit succeeds 
in reproducing the fact that the data at very small t are higher at low energies than at the Tevatron. 
Figure 5b shows our expectations for LHC energy; these predictions should be regarded as qualitative 
because, as we have explained, there is not a sufficient quantity of good lower-energy data to make a 
definitive fit. 

Figure 6a shows da/dt. Although our fit used only data for M 2 > 9 GeV 2 , we have extrapolated it 
down to 1.4 GeV 2 and integrated over £ down to this value and up to £ = 0.05. We have multiplied 
by 2 to account for both pp — > pX and pp — ► Xp (or pp — ► pX and pp — > Xp). We show the plot only 
down as far as —t = 0.036 GeV 2 , for the reasons we have explained. We have explained also that it is 
risky to extrapolate the simple fit to values of M below M = 3 GeV. 

Figure 6b shows the integral of da/dt. The middle curve integrates it down to —t = 0.036 GeV 2 and 
the lower curve shows how much comes from the term 1P p . We have extrapolated the curves in figure 
6a down to t = to get a rough estimate of how much we should add to the middle curve to make 
it extend to t = 0. This extrapolation is highly uncertain. If we assume that between t = and 
—t = 0.054 GeV 2 da/dt may be well approximated by a single exponential e bt , so that 

<r(t M1N = 0) = (<r(t MIN = 0.036)) 3 /HtM IN = 0.054)) 2 (18) 

this results in the upper curve. If instead we had used the parametrisation (17) with our large negative 
^rrr it would have led to a divergent cross section. 

The upper curve in figure 6b fits well to 

<t dipf (s; M 2 IN = 1.4 GeV 2 ) = 1.9 s ' 1425 (19) 

in mb-GeV units. The data are from reference [15]. We must emphasise that experiments all have 
limited acceptance in t and £, so that the values they quote for <r DIFF (s; M MIN ) are more than a little 
model-dependent. As we have just explained, our fit to the available data, which are for —t > 0.036 
GeV 2 , can give very different total "diffractive" cross sections when it is extrapolated to t = in 
different ways. Because we do fit the available differential-cross-section data, the total diffractive 
cross section represented by the upper curve in figure 6b is just as believable as the claimed data 
points; there is no trustworthy extraction of the total diffractive cross section. This in turn leads to 
uncertainty in the results for cr TOTAL when it is extracted from collider experiments by the so-called 
luminosity-independent method t 2l l . 

Our fit (6) to the total cross section predicts cr TOTAL = 101.5 mb at yfs = 14 TeV, while (19) gives 
^diff _ 29 m j-^ so t nere i s no need to introduce significant shadowing up to this energy. 

As is seen in figure 6b, even at Tevatron energy the term F p contributes less than 70% of the 
total "diffractive" cross section. At very large energies this fraction rises to nearly 80%. Accurate 
measurement of d 2 a/dtd^ at RHIC would be of great value in clarifying further the details of the 
diffractive dissociation mechanism. 
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